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We give a detailed renormalization group analysis for the SU(16) grandunified 
group with general breaking chains in which quarks and leptons transform separately 
at intermediate energies. Our analysis includes the effects of Higgs bosons. We show 
that the grandunification scale could be as low as ~ 10^'^ GeV and give examples 
where new physics could exist at relatively low energy 250 GeV). We consider 
proton decay in this model and show that it is consistent with a low grandunification 
scale. We also discuss the possible generation of a neutrino magnetic moment in the 
range of 10^^^ to 10^^^ fiB with a very small mass by the breaking of the embedded 
SU(2),y symmetry at a low energy. 
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I. INTRODUCTION 



Recent treatments of the SU(15) grandunification group [jlj-Q have shown that in the ver- 
sions of that model which produce "ununified models" at intermediate energies, granduni- 
fication may be reached at a relatively low energy and that the lowest intermediate scale may 
be within the reach of the Superconducting Super Collider (SSC). It has also been shown 
recently that the effects of Higgs bosons in the renormalization group equations of SU(15) 
models are quite large [Q. In the literature, there is also some discussion of the gauge group 
SU(16) 10], which has some desirable features that are not found in SU(15) models. Among 
these are seperate gauging of baryon and lepton number and the embedding of Voloshin 
symmetry P,^] Sl]{2),y which might play an important role in the solar neutrino puzzle 
||TO| , p]T|| . In this paper, we will analyze the renormalization group equations (RGE's) of the 
SU(16) grandunification model with breaking chains in which quarks and leptons transform 
separately at intermediate energies. We will include the effects of Higgs bosons, which are 
significant, and use the newest LEP data for couplings at low energy. We give this analysis 
in Section 0. 

In the SU(16) model, all known left-handed fermions of a single generation together with 
a left-handed antineutrino transform like the fundamental representation of the gauge group. 



I.e. 



= (^'e i^e e 6^ Ui U2 Uj, di 62 dj, Ul U2 di ^ (1.1) 

Here, the hat symbol over a particle's symbol denotes its antiparticle. Of course, mirror 
fermions must be introduced to make the model free of anomalies, but we do not need to 
discuss them explicitly here. Our interest is to look for chains with low unification scale. 
Existence of such chains is known [|]-^ in SU(15) models, for which it has been shown that 
a low unification scale makes the model free from the monopole problem [|1^ while being 
perfectly consistent with the proton lifetime |p|,pr3|. We discuss proton decay for our SU(16) 



model in Section |IT1|. 
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Another reason for examining the SU(16) model is to implement the suggestion of a 
previous paper which points out that this group contains the subgroup SU(4)/, in which 
the left-handed leptons including an antineutrino transform as the fundamental representa- 
tion. Further, 811(4); contains Voloshin symmetry which allows for the magnetic moment 
of the neutrino to be in the range of 10^^^ to IO^^^^b and at the same time allowing only 
a small neutrino mass P, [l4| , [T5| . Magnetic moment in this range might be required |]16| if 



the anticorrelation |jTT| of solar neutrino flux with sun spot number is confirmed. We shall 



discuss the implementation in detail in Section |y 



II. SU(16) BREAKING SCHEME AND RG ANALYSIS 

In Fig. we show the symmetry breaking scheme of our model. The purpose of our 
scheme is to have the leptonic sector transform separately from the quark sector at inter- 
mediate energies and to have the leptonic sector SU(4)i break in the same manner as in the 
previous SU(4); based model 0. In the figure, we show the representations of the Higgs 
fields used to break symmetries at their indicated mass scales. We use the hypothesis of 
minimal fine-tuning , which allows us to choose the mass scales of submultiplets of these 



fields at different scales. We are interested in scenarios where an intermediate gauge group 
exists at energies of the order of 1 TeV so that we shall have new physics at observable en- 
ergies. We will show in this section that such scenarios are consistent with renormalization 
group analysis. 

In the one-loop approximation, couplings ctj = gf/4:n evolve as 

dai Bi 2 

-7^(^i , (2-1) 



dlnM 271 
which implies 

a;\M2) = ar\M,)-^\n^^. (2.2) 

For the standard model couplings aiy, a2L, and a^c, we use the conventional normalization 
which determines the relations 



Tr(T,T,) = 25,,, = - ^{2ng) - ^T{S,) , (2.3) 

where the Tj's are the generators of the SU(16) fundamental representation. To simphfy the 
boundary conditions, we normahze the couphngs of the intermediate gauge groups and 
the SU(16) gauge group such that 

Tr (T,T,) = , B, = j (^^N - ^(2n,) - ^T{S,)^ . (2.4) 

Here, / is the number of fundamentals of the subgroup per generation in the 16 of Eq. ( |1.1D . 
To be explicit, for the SU{2)gi group / = 3, for the SU(3)qL group / = 2, and for all other 
intermediate groups and the SU(16) group / = 1. The above equations for Bi hold for 
SU(N) gauge groups and with = hold for U(l) gauge groups. In the second terms, Ug, 
the number of fermion generations, is multiplied with a factor of two to account for mirror 
fermions. This term does not affect the scales of symmetry breaking. In the third terms, 
T{Si) is the quadradic Casimir invariant for all Higgs submultiplets with masses less than 
the scale of interest. For a complex field, the value of T{Si) should be doubled. 

With the above normalizations, the U(l) generators that enter our symmetry breaking 
pattern (see Fig. |l|) follows: 



QqB = 




(0(4), 1(6), -1(6)) , 


(2.5) 


QqA = 




(0(4), 0(6), 1(3), -1(3)) , 


(2.6) 


QqY = 


— 1= dia; 
2V33 


?(0(4), 1(6), -4(3), 2(3)) , 


(2.7) 


Qix = 




(1,1,1,-3,0(12)), 


(2.8) 


QlY = 


2V3''"'* 


(0,-1,-1,2,0(12)), 


(2.9) 


Qy = 


— = dia; 
2^15 


g(0, -3, -3, 6, 1(6), -4(3), 2(3)) . 


(2.10) 



In these equations, notation of the form a(b) stands for b successive entries of a. For example, 
"1(6)" stands for "1,1,1,1,1,1." The ordering is the same as in Eq. (|1.1|). 

The breaking scheme shown in Fig. |I] has the following boundary conditions at the 
breaking scales Mf. 
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— 1 / 71 J- \ 


«ri(^i2) = "6i(^^i2) 


= (^gMMu) = 




— 1 / 71 /T \ 






1 1 

^^3ur{Mb) + ^asdRiMB) 




2 q 

-—1 / 71 ,r \ 1 — -1 / 71 ,r \ 

— aig^B(MB) + — «i^(Mb) 


— 1 / 71 ,r \ 


— 1 / 71 /T \ 


— 1 / 71 /T \ 

= (^31 (Mil) = c 


— 1 / 71 /T \ 


— 1 / 71 /T \ 

= «2iL(^3i) , 


1 ^ 

— 1/71^^ 1 — l/Tl/TN 

gttsi (^3/) + -ai/x(^3i) 


= (^iiviMsi) ; 


2"3LV^^y) + 2«3i(^y) 


1 / 71 J- \ 

= "3c (^y) , 


3a2"i(My) + a2KMy) 









hqBi 



-1 



(2.11) 
(2.12) 
(2.13) 
(2.14) 

(2.15) 
(2.16) 

(2.17) 



(2.18) 



In writing Eq. ( |2.17| ), we have used the fact that in the group SU(3); there is a generator 

1 



A 



31 



2^3 



diag (2, -1,-1, 0,0(12)) 



(2.19) 



which is broken at the scale M3;. At the same scale U{l)ix also breaks, leaving unbroken 
the combination 



QiY — -A3/ '^Qix 



(2.20) 



In deriving Eq. ( |2.18|) for the breakings at the scale My, we have used similar considerations. 
In particular, SU(3)qL, which has the two diagonal generators 

1 



3qL = diag (0(4), 1, -1, 0, 1, -1, 0, 0(6)) , 

^ diag (0(4), -1,-1, 2, -1,-1, 2, 0(6)) 



breaks, as does SU(3)gR, which has the two diagonal generators 



(2.21) 
(2.22) 
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Asgii = ^ diag (0(4), 0(6), -1, 1, 0, -1, 1, 0) , (2.23) 

^3qR = ^ diag (0(4), 0(6), 1, 1, -2, 1, 1, -2) , (2.24) 

leaving unbroken SU(3)c which has the two diagonal generators 

Aac = V2X3qL + V2X3qR , (2.25) 

A3, = V2A3qL + V2A3gi? ■ (2.26) 

Similarly, SU (2)^^ with the diagonal generator 

^2qL = diag (0(4), 1(3), -1(3), Oe), (2.27) 
and SU{2)iL with the diagonal generator 

A2iL = ^ diag (0,1, -1,0, 0(12)) (2.28) 
are broken, leaving SU{2)l unbroken with the diagonal generator 

A2L = V^A2,L + A2/L- (2.29) 

Now, using 

^ logio (t^) , (2-30) 



''HlGeV, 

the one-loop equations and boundary conditions give us the following equations for the 
standard model couplings: 

ZTT 

+ (ri4/ - n3z)(2S3,i + 2^3,^) + {ub - n4l)i2B3gL + 2B3gR) 

+ {nQR - nB){2B3gL + B-^uR + B^r) 

+ {n6L - ri6il)(2-B3qL + 2^6^^) 

+ (ni2 - neL)i2BegL + '^Bq^r) 

+{nG - ni2)4Bi2q] , (2.31) 



a^liMz) = Aaa^ - ^^iiny - nz)B2L + - ^y)(352gL + B21L) 
Zn 

+ {n4i - n3i){3B2qL + B31) + {hb - n4i){3B2qL + B41) 

+ {n6R - nB)iW2qL + B4i) + {uql - neR){3B2qL + B^i) 
+ (ni2 - neL){^BQqL + B^) 

+ inG-nu)iWi2qL + B4i)], (2.32) 

aj"y(Mz) = 4a5^ - -^[(it-y - nz)BiY + (^3; - r2y)(y5igy + '^Buy) 

11 1 8 

+ {nii - n^i){—BiqY + -^3, + -Biix) 
5 5 5 

11 9 

+ {nB - n4i){—BiqY + -B41) 

2 9 9 

+ (^6fl - nB){-BiqB + g^lAg + ^^40 

2 9 9 

+ {nQL - nQR){-BiqB + -^BQqR + -B^i) 



z y y 

+ (^12 - nQL){-BiqB + g^egi? + g54«) 

+ {nG - ni2)(ySi2, + ^54z)] . (2.33) 



For the couphngs at Mz we use the experimental values |TB 



a3-^(Mz) = 8.197, 
a2-2(A^z) = 30.102, 
a-^(Mz) = 59.217, 

= 91.176 GeV. (2.34) 



The BiS are determined by Eq. ( |2.3| ) and Eq. ( |2.4| ) with the T(S'i)'s being determined by 
the Higgs structure given in Fig. |l| along with the principal of minimal fine-tuning. We note 
that to impliment the suggestion of ref. an additional rank-2 antisymmetric tensor is to 
be included. However, this Higgs field has little effect on the RGE's and is not included in 
the analysis of this section. 

Now, the above equations for the standard model couplings can be solved silmutaneously 
in terms of hq and ni2 to obtain 



no = -3.28 - 0.09ny - 0.21^3/ + Omn^i 

+0.32nB + 2.99nQR - 2.60n6L , (2.35) 
ni2 = -2.08 - 0.02ny - 0.22n3, + O.rSnn 

+0.37nB + 2.70nm - 2MnQL ■ (2.36) 

We note that Higgs fields make significant contributions to the above equations. From the 
structure of the above equations, we see that n^i and uqr being relatively high helps to meet 
the requirement no > ni2. In fact, no solution with low 71,4; is found to be acceptable. 

We are further interested in seeing that it is possible to have new physics, including 
breaking the Voloshin symmetry, at less than 1 TeV. Since our lowest intermediate stage, 
which contains the Voloshin symmetry, is broken at M3;, we are interested in n^i ^ 3. So as 
to investigate this possibility, we make the simplification 

M^L = M12 = Mg , 

M^i = Mb = Men . (2.37) 

This yields 

nc = 9.35 + 0.66ny - 0.28n3i , (2.38) 
ngR = 8.77 + 0.59ny - 0.19n3i . (2.39) 

In this example both My and M^i may be low. The particular case of this solution which 
has My = M3; is graphed in Fig. ^ for the allowed region where Mq > M12 > M31 ^ 250 
GeV. This gives the range of unification scale to be 10^° GeV < Mg < 10^^ GeV. 
We also investigate another possible solution where 



Mil = Mb = Men = Mql = Mu 



My = 10^ GeV . (2.40) 

This yields 
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na = 3.99 + O.SGna/ , (2.41) 
ni2 = 4.62 + 0.48n3/ . (2.42) 

We graph this case in Fig. ^. Note that Mq can be as low as 10^'^ GeV. In fact, from the 
constraint M12 > M3/, we note from Fig. ^ that the unification scale has to be smaller than 
10^'^ GeV. One characteristic of this solution is that all intermediate scales other than My 
are larger than 10^ GeV, so that the only observable new physics comes from the scale My- 
However, we shall see in the next section that this solution is inconsistent with the constraints 
arising from nonobservation of proton decay, unless the discrete symmetry —V^ is 

imposed on the Lagrangian. 

Each of the solutions discussed above has the feature that one can 
achieve grand unification at scales much lower than what is expected in standard uni- 
fication models based on gauge groups SU(5) or SO(IO). Such low grandunification scale 
has many interesting features not present in standard unification models. For example, it 
has been shown |jl2| that, unlike the SU(5) model, one can circumvent the cosmological 



monopole problem easily. Also, the intermediate scales are low, and some of them can be in 
reach of the next generation of experiments. If new gauge structure exists at TeV energies, 
as occurs for example in the second solution with low My, there will be new gauge bosons 
to mediate a lot of processes 0. Direct production of these gauge bosons at the SSC will 
also provide exciting physics, some of which has been discussed in the context of SU(15) 
models [0,0. 



III. PROTON DECAY 

Since baryon and lepton numbers are part of the gauge symmetry of the model, proton 
decay diagrams must involve baryon number and lepton number violating vacuum expecta- 
tion values (VEV's) in the Higgs sector. The only VEV in our symmetry breaking scheme 
which violates baryon number is the VEV of the 560-dimensional representation, B^^"^, hav- 
ing the quantum numbers of (uddj, which has B = —1. Lepton number violation arises 
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from the VEV of i?'^'™ having the quantum numbers of {d{ue — dve)ji which has L = 1, 
and the VEV of 16-dimensional representation, V^", having the quantum numbers of (z>e), 
which has L = —1. Note that lowering indices changes the signs of the quantum numbers. 

To examine proton decay, we use the method of effective operators. In this method, one 
writes down all effective operators involving fermions and scalars which are invariant under 
the full SU(16) group. When the scalars develop VEV's, one obtains an effective operator 
involving fermions only. These VEV's are responsible for baryon or lepton violation. 

The lowest order effective operators for proton decay involve four fermionic fields P|. 
From the discussion above, we find two SU(16)-invariant effective operators which can induce 
proton decay |2^. These are 

= {^'=^'}{^"'^"}5fc,,5P'^^$„p$„,, , (3.1) 

O2 = {^"^ji^^'^^jBklmVn . (3.2) 



Here, {^*^^'} = y^'^j C^i, where C is the conjugation matrix for fermions. 

Typical diagrams generating Oi and O2 are shown in Fig. ^ and Fig. ^, respectively. 
In Fig. ^, ^uii gets a VEV. In determining the RGE's, we only gave a VEV to the $^ 
component of 136. However, in order to give masses to the quarks as well as the charged 
leptons a VEV must be given to each of the and $^ '^^ components of 136, i.e. a 

linear combination of and represents the standard model Higgs. At higher 

scales additional components make contributions to the RGE's. These contributions are 
small compared to other Higgs contributions. Therefore, we ignore them. 

From the figures we obtain an order-of-magnitude estimate of the coeficients ki and K2 
of the 4-fermion operators Oi and O2 respectively: 

mjY Xb^XbbMyMbM^ 



\MwJ ' ^ ^ 

Here, the quantity is the mass of a typical fermion, and comes from the Yukawa couplings. 
Antisymmetry of B^ir and fermion indices require use of second generation fermions [^. 
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Therefore, we use mj ^ 100 MeV. Also, Xb<s>, ^bb, and A denote the scalar couplings. We 
have assumed that all virtual colored scalars have masses of order Mq, the largest scale of 
this model. The mass scales in the numerator are the scales of the VEV's. Making a rough 
estimate, we have neglected factors of gauge coupling constants. 
Known bounds on proton lifetime imply 

Ki , K2 < GeV-\ (3.5) 

If we use the above constraint with Mb ~ Mq, which can be seen to be true from Section 
H, and assume Ab*, Xbb, and A are ~ 1, then we find the constraints 

^ > 10^^ GeV^ , (3.6) 

^>10'^GeV'. (3.7) 
M31 

from Oi and O2 respectively. Eq. ( p.6|) puts no restriction on the solutions to the RGB's, 
but Eq. ( |3.7|) , resulting from O2, does. For example, Eq. (|3^ ) rules out the entire special 



case defined by Eq. (|2.40|) which gives Mq as low as ~ 10*^'^ GeV, although it does not rule 
out any region of the case defined by Eq. (|2.37|) which allows for a low energy M^i. However, 



the effective operator O2, which yields this constraint, is no longer allowed if we impose the 
discrete symmetry V" —V" on the Lagrangian. Another feature that would exist if this 
discrete symmetry is imposed is mentioned in the next section. 

It is important to note that the decay modes of the proton obtained from the operators 
Oi and O2 are different from the predictions of standard SU(5) or SO(IO) grandunification 
models. In Oi, the indices k,l are antisymmetrized and so are m,n. Thus, the quark level 
operator for proton decay ^ arising from it is usufi~^. This gives rise to the decay mode 
p — > K^fi~^. On the other hand, in O2, the quark level operator is usdi), where the neutrino 
can belong to any generation since the indices m, n are not necessarily antisymmetric for 
this operator. Thus, we expect a decay mode p K^v. Note that both operators give rise 
to B — L conserving decays. 
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IV. NEUTRINO MAGNETIC MOMENT 



We now discuss generation of a sizeable magnetic moment for the neutrino. The most 
general Yukawa couplings of the model are given by 

-Cy = Y. haKA^cP + E fabKA^cP + h.C. (4.1) 

Here, Latin indices refer to the generation. $ is the symmetric 136-dimensional rank-2 
SU(16) tensor representation, whose couplings in the generation space are chosen diagonal 
without loss of generality. The additional multiplet (f is a. 120-dimensional antisymmetric 
rank-2 tensor, so its coupling fab is antisymmetric in its generation indices. This field is put 
into the model to generate a magnetic moment for the neutrino. Since quarks play no role 
in the magnetic moment of the neutrino, we focus our attention on the leptonic part of the 
interactions. 

$ is the Higgs field which breaks SU(3)c x SU(2)^ x U(l)^ down to SU(3)c x U(1)q. In 
the leptonic sector we assume only $34 gets a VEV 0. This VEV gives masses to charged 
leptons, but not to neutrinos. The multiplet (f, on the other hand, is assumed to have 
no VEV [^]. The dominant contributions to the mass of z/g then come from the one- loop 
graphs of Fig. |^. The diagrams involve a Higgs potential term of the form '-y{p°'^{p'^^Aa/3:ys, 
where 7 is a coupling with the dimension of mass and A^i^^s is the antisymmetric rank-4 
SU(16) tensor whose VEV (A1234) breaks SU(16) to Sf/(12), x SU{A)i at Mg- Because 
A-ap^s is antisymmetric, (A2314) = — (A2413), and therefore the mass contributions of the 
two diagrams of Fig. H have opposite sign. 

From Fig. |^, we estimate 

where Mq^ = (massof y^a/j). The important point to realize here is that, in the limit of 
unbroken SU(3);, M23 = M13 and M24 = M14, and so the mass contributions from the two 
diagrams cancel each other. 
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The diagrams of Fig. ^ with one photon hne attached give the dominant diagrams for 
the magnetic moment of v^.. These diagrams add because of an extra negative sign between 
the two diagrams arising from the photon vertex. We estimate 



(^fisl (^2413) rrir 
167r2 



1 



1 



(4.3) 



In the hmit of unbroken SU(3)i, we can write M'^ 
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MI3 = M2 - \AMl and Mf^ 



MI4 = M^ + ^AMli, where AM|; is the mass difference due to SU(4)/ breaking. The breaking 
of SU(3)i then changes the masses according to 



- ^AM^ + AMI , Mi, = M'- ^AM^ , 
Ml = M' + ^AMii + AMI , Ml, = M' + ^AM^ , 



(4.4) 
(4.5) 



where AM|; to lowest order is due to the term A (Vi) {V ) (f (pia- Here, Va is the vector 



representation Higgs field and A is a dimensionless coupling. Putting Eq. ( |4.4| ) and Eq. 



( [4.5|) into the expressions for m^^ and /ij,^ gives us 



f^slMgm 
167r2M6 3« 



41 



'AMlX 



2M2 



effslMgnir 
Stt^M^ 



1 - 



' AMI/ 
2M2 



(4.6) 
(4.7) 



where we have used (A2413) ~ Mq and assumed <C 1. Now, requiring m,^^ < lOeV and 
fJ'Ue ~ lO^^^yUs, gives from Eq. ( ^4.6| ) and Eq. ( |4.7| ) the constraints 



/f37M, 



^AM|,AM|, 
M4 



'am: 



2 \ 2 



4/ 



2M2 



1 



'AM 



4/ 



2M2 



< 10"^ 



> 10^6 GeV-^ 



(4.8) 
(4.9) 



Also, demanding the mass-mixing matrix for (f to have a positive determinant gives the 
constraint 



7 < 



M^ 



M 



G 



1 



' AMI/ 
2M2 



(4.10) 
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We remark that if we impose the discrete symmetry — — on the Lagrangian, then we 
should expect 7 to be small because then '-yip°'l^ip"'^ A^js-ys is the only term of the Lagrangian 
not invariant under —Aap-yS- An example of this is when in each of Eq. ([4.8D, 

Eq. ( |4.9| ) and Eq. ( |4.1CI| ) the left and right sides are approximately equal. In this case, 
AMliAMii ~ lOVfsGeV^ and if also (^^^ < 1 then M ~ lO^ GeV. Typically, M is 
of the order of ITeV with SU(3)/ and 811(4)^ breaking inducing small mass changes in the 
ipap's with AMil ~ AMI ~ 10=^ GeV^ 



V. SUMMARY 

We have shown in this paper breaking chains for the SU(16) grandunification group which 
lead to the standard model with the quark and lepton sectors transforming separately at 
intermediate energies. The grandunification scale for this model can be as low as ~ 10^'^ GeV. 
We have shown that this does not produce any conflict with know bounds of proton lifetime. 
In fact, if a discrete symmetry is imposed on the model, one can obtain chains with the 
unification scale as low as ~ lO*^'^ GeV. Also, low intermediate breaking scales can exist in 
the ^ 1 TeV range. This has many observable consequences, including guage bosons at the 
TeV range which can give rise to a rich phenomenology. 

Further, the model embeds the Voloshin symmetry SU{2)i^ into its subgroup SU(4)/. By 
using a rank-2 antisymmetric Higgs field, it is possible to get a significant magnetic moment 
for a neutrino with small mass. For this to work, this Higgs field typically should have masses 
~ 1 TeV with relatively small mass differences induced by 811(3); and 811(4); breaking. 

This work has been supported by the Department of Energy grant DE-FG06-85ER-40224. 
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FIGURES 

FIG. 1. A possible chain of symmetry breaking. The numbers n denote a factor SU(n) in the 
gauge group if n > 1, and a U(l) factor if n = 1. The superscripts q or / indicate whether only 
quarks (and antiquarks) or only leptons (and antileptons) are non-singlets under that part of the 
gauge group. If one considers the 255 as a traceless matrix, its VEVs are diagonal and the notation 
1(g), e.g., stands for six consecutive entries of unity. In the 560, the symbol (widehatdue) , e.g., 
stands for the VEV of the color singlet combination of the components with one index having the 
quantum numbers of d, another of u and another of e. One can contemplate chains with fewer 
steps by equating two or more energy scales. Note the following transformation properties: A^^^''''^ : 
1820, TI' : 255, i^f^^j : 14144, B^'^''^ : 560, V : 16, <^>^'^^ : 136. 

FIG. 2. no = logio w ^^'^ "-6/? = logio "Sif ^ ^ function of n-M = log^o with My = M31, 
M^i = Mb = Mqr and Mql = -^12 = in Fig.l. The shaded area is not acceptable since by 
definition Mq > Mqr > M31. 



FIG. 3. no = logio and nu = logio ^ ^ function of 713/ = logio "^^^^ 
Mil = Mb = Mqr = Mql = M12 and My = 10^ GeV in Fig. 1. The shaded area is not 
acceptable since by definition Mq > M12 > M^i. 



FIG. 4. Tree diagram generating Oi. The labels on the Higgs boson lines represent, via Eq. 



(1.1), the transformation properties under SU(16). The indices should all be considered as upper 



indices. 



FIG. 5. Tree diagram generating ©2. 

FIG. 6. The diagrams that contribute to the mass of the neutrinos at the 1-loop level. Magnetic 
moment arises by attaching a photon line to any internal line. 



17 



